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ABSTRACT 


The  analytical  treatment  of  the  problem  of  solidification  in  a slab 
of  thickness  2a  is  based  on  Lightfoot's  solution  for  the  case  of  a semi- 
infinite  mass,  initially  at  constant  temperature  -vnd  bounded  by  the  plane 
surface  x — 0 kept  at  constant  temperature.  The  rate  of  advance  of  the 
solid  walls  moving  from  the  two  boundary  planes  x — 0 and  x = 2a,  kept 
at  constant  temperature,  is  given  in  closed  form.  Solidification  in  a 
sphere,  in  an  infinite  or  finite  cylinder,  or  in  a parallelepiped,  is 
obtained  from  the  solution  for  the  slab  by  correlation  of  temperatures  due 
to  the  initial  supply  of  heat.  Numerical  solutions  can  be  calculated 
rapidly. 


1.  INTRODUCTION: 


An  analytical  solution  of  the  problem  of  freezing  in  a finite  region 
has  been  attempted  by  Lightfoot#  who  studied  the  rate  of  solidification  of 
a medium,  originally  at  constant  temperature,  and  bounded  by  the  plane 
surfaces  x = 0 and  x = 2a  kept  at  constant  temperature.  The  problem 
was  treated  by  the  method  of  images  but  essential  factors  were  neglected 
and  the  solution  obtained  by  Lightfoot  is  neither  complete  nor  exact. 

An  exact  solution  has  been  obtained  by  the  same  author  for  the  case 
of  a semi-infinite  mass  of  molten  material,  originally  at  constant  tempera- 
ture, and  with  its  bounding  plane  surface  x,  = 0 maintained  at  constant 
(zero)  temperature.  In  the  given  solution,  the  diffusivity  k,  the  specific 
heat  c,  and  the  density  p of  the  medium,  are  assumed  constant  for  all  tempera- 
tures, and  the  same  for  the  liquid  as  for  the  solid  medium. 

The  so'ution  of  the  problem  of  solidification  in  finite  regions, 
given  below,  has  oeen  derived  from  Light  foot's  solution  for  a semi-infinite 
medium.  The  following  method  will  give  the  position  of  the  surface  of 
separation  of  solid  and  liquid  phases  as  a function  of  tine.  Results  of 
this  study  may  be  of  help  in  the  investigation  of  a complete  so'ution  of  the 
problem  of  solidification. 

2.  Solidification  in  a Semi-infinite  Medium**- 

The  solution  of  the  problem  of  solidification  of  a mass  of  molten 
material  requires,  in  addition  to  known  solutions,  consideration  of  the  effect 
on  the  temperature  caused  by  the  evolution  of  latent  heat  of  fusion.  In 
Lightfoot's  analysis  of  the  problem,  the  surface  of  separation  of  solid  and 
liquid  phases  has  been  regarded  as  a moving  source  of  heat.  The  speed  of 
of  the  plane  source  was  determined  from  the  condition  that  the  temperature 
at  the  moving  plane  surface,  at  any  position,  is  equal  to  the  melting  point 
of  the  material. 

The  position  x of  the  moving  plane,  at  time  t,  in  a medium  initially 
at  temperature  and  with  its  bounding  plane  x — 0 kept  at  zero 
temperature,  is  given  by 

x = 2*.</kt  (1) 

where  k is  the  thermal  diffusivity  of  the  medium,  and  the  constant  jC  is 
obtained  from 

v • p e (*)  + HS*  t<e^*e  (*)  (i  -©(*)]  (2) 

* Proc.  London  Math.  Soc.  (2),  31,  (1930),  97 
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terms  used  in  (2)  are  as  follows: 
the  melting  point  of  the  medium 
the  initial  temperature 
the  latent  heat  of  fusion 
the  specific  heat 


(3) 


With  this  solution,  the  temperature  V,  at  a point  xt,  at  time  t, 
for  x,  <z4t/JEI  will  be: 


f**e(zvfeF)+  (z^t)  {'-«(*)} 

and  for 


According  to  (4)  ar.d  (5),  we  notice  that,  while  solidification 
proceeds  from  the  plane  x — 0,  the  temperatures  due  to  the  original 
sunply  of  heat,  and  to  the  latent  heat  of  fusion,  can  be  evaluated  separately. 

For  x = at  the  surface  of  separation  of  solid  and  liquid 

phases,  (4)  and  (5)  coincide  with  (2). 

The  terms  in  (2)  are  constant  for  any  position  of  the  moving  plane 
surface  source,  and  can  be  written  as: 


V,  s *©(*) 

VL  • ^ <«**©(*)  (•-©(*)} 

V - vf  + vt 


(6) 

(7) 

(8) 


3 


\ 


3.  Solidification  in  a Slab  of  Thickness  2a 


An  exact  solution  of  the  problem  of  solidifies  ■ : " ' ' i »•  a*  of  molten 
material,  initially  at  constant  temperature,  can  ■>>.  derived  from  the  results 
given  in  the  preceding  section. 

For  a freezing  medium  in  the  region  2a,  initially  at  temperature  $ , 
with  the  bounding  plane  surfaces  x = 0 and  x a 2a  kept  at  zero, 
we  evaluate  the  sum  of  the  temperatures  due  to  the  original  supply  of  heat 
and  to  latent  heat  of  fusion  released  on  solidification  by  the  two  solid 
sections  of  the  medium  (0,  x)  and  (2a,  2a  - x)  as  in  Fig. 1. 


Fig.  1 

The  surfaces  of  separation  of  solid  and  liquid  phases  are  taken  as 
two  moving  ’■'lane  sources  of  heat,  and  the  speed  of  these  two  sources  is 
determined  from  the  condition  that  the  temperature  at  the  moving  plane 
sources  is  equal  to  the  melting  point  of  the  medium. 

For  too  tennaralrro  at  the  moving  sources  we  write  the  condition 

V=  v^  + Cv-  )j  + (v)2 

where  refers  to  the  original  supply  of  heat,  and  (Vj_  , 

are  temperatures  due  to  heat  of  fusion  from  solid  section  (0,x) 

(2a.  2a  - x)  respectively. 


(9) 

(VL  )Z 
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It  is  evident  that  solidification  will  proceei  from  both  x - 0 and 
x n 2a  , and  for  t snail,  the  positions  of  the  moving  planes  will  be 


X - 2 f 2* -X  = ZA-Zt'VXTt 


(10) 


The  constant  ^ of  (10)  is  the  same  as  for  a semi-infinite  medium  and 
is  given  by  (2).  In  terns  of  the  dimensionless  variable 


forms  (10)  will  become 


£*r.JL-  z~  2*#-/r 

s 


(li) 


(12) 


In  the  region  the  temperature 
given  T"  , is  the  known  solution 


in  + l 


s in 


(Zn  + l)  nx. 

2T3 


, at  any  point  x , at  a 

-(Zn+l)Zn* 

6 


(13) 


During  an  initial  period  of  time,  (V  ) and  (V^)  will  be  obtained 
from  ( O and  (5)  as  follows:  * 


For  |t<24#vT‘(VL 

>,■ 

LV*rc 

e 

6 

(lSvw)\ 

'1-  8(0} 

(14) 

and  ( Vj_ 

L-ffi 

c 

* * 

1 

© 

'am) 

(15) 

For|t>2^i/T(VL 

1 = 
' 1 

L-fiE 

c 

X * 

A*  e*#  e 

60  {'-0 

6dhr)} 

(16) 

(vL 

same  as  (15) 

For  the  moving 
be  evaluated  from: 

plane  source  — - s ZA0  yx  t 

.he  terms  in 

(9)  will 
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_ . . »irv  _ 

ZH+ I 24 


v>  s ) ' ! 

f n n?c 

(VL\ » *=?  <.«*•*©  «)  (|-  e (-a.)J 

(vL)t- HF  *•***  ®«>  I1- 9 (lIvT-)) 


I (2n+lWlx  -(2n  + l)  TT  t4 

y , ®in  — c 


(17) 

(18) 

(19) 


^ Tempe rature  , (V^  ) , and  (V^Ljat  the  moving  plane  source 

will  be  obtained  with  U-7),  (19),  and  (18),  respectively. 
The  sum  of  (17),  (18),  and  (19),  is  equal  to  the  melting  point  of  the 
medium  for  a range  of  values  O 4 T 4 . The  upper  limit  7^  varies 

for  each  numerical  case  and  depends  upon  the  initial  temperature  p and 
the  ratio  -j?-  . For  7"  srn^ll  , the  temperature  at  the  center  of 

the  slab  is  given  by 

v-  _ 4^  I (Zn+l)Tl  -(Zn+l)*'*TZ*‘ t/4. 

x-4  -=rL  7TTT-7  Sln  „ / e 


+ 2 ^ V4**9  «.>  { I - • (*7?)} 


(20) 


and  we  can  verify  that,  when 

T^T-  »<*-»)*  ^ 


In  each  numerical  cas  ),  therefore,  the  upper  limit  7^  can  be  obtained 
with  {20).  We  see  also  that,  during  the  entire  period  of  time  *U9  , 

the  boundary  conditions  of  the  problem  will  be  satisfied. 

..  In  the  case  of  a semi-infinite  medium  and  x > 2*V^F  , or 
■3T  >2^0  VT  , the  temperature  given  by  (5)  is  the  sum  of  the 
two  terms 


(21) 

(22) 
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Solution  (21)  is  correct  for  large  values  of  *T  but,  for  small  values, 
must  be  equal  to  f at  a finite  distance  frar.  the  boundary 
x = 0.  Where  , we  will  obtain  (^i,)9  0.  In  general,  the 

boundary  condition  is  satisfied  when  the  temperature  at  the  center  of  a 
slab  of  thickness  2* a , initially  at  0 and  with  bounding  surfaces  kept 
at  zero  is 


I 


2 e (23) 


<VAc  . 24 -*  = Sm 

t«t?  n*c 

This  condition  can  be  verified  for  any  70  , including  the  limiting  case 
where 


X*  a 2 ipt  » I 

* 


and  we  can  write 

(VLW“° 

r 

During  the  period  of  time  L0  , the  two  plane  sources  move  from  the 
boundaries  toward  the  center  of  the  slab  as  in  semi-infinite  medium. 

The  temperature  distributions  , (Vj_  ),  and  (V^_  L and  their  sum, 

computed  for  a numerical  example,  with  T*  , are  shoWn  on  Fig,  2. 


In  the  region  2a,  initially  and  till  solidification  is  complete, 
the  correct  mathematical  forms  of  (Vj_  )f  and  (*u  , at  the  plane 
source  moving  from  x = 0 to  x = are: 


n«l  J0 


sin  n1tx'  'n‘n  hx'.T>(2h) 

sm~z T°  3fdt 

T-T'l 


n»i 


2« 


i n7rfea»-x')  "n*7r  ^ 4 ^ ■> 

"7  s z3 — * S?dT  (25) 


During  the  interval  o<tT<76,(vL)  is  a constant  obtained  from  (18) 
while  (V  j_  )fc  , evaluated  from  (19)  , is  initially  zero  but  may 
assume  increasing  values  as  T approaches  T0  . However,  during  the 
entire  interval  O £ Tit 

(VL)t  . v-{v,  + (vu),}  (26) 
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The  rate  of  solidification  in  the  range OjgJC^^mi ght  be  obtained 
directly  from  (24)  and  (25)  if  the  series  of  integrals  in  the  above 
forms  oould  be  evaluated  exactly.  So  far,  available  graphical  or 
numerical  techniques  are  unsatisfactory. 


Light  foot  shows  that  a rate  of  solidification  can  be  obtained  by 
trial  method,  and  assumes  a position  of  the  moving  planes  given  by 


b-2Ky^f 


T>TC 
*2  *© 


where  b,  K,  and  , are  chosen  for  continuity  and  satisfy  the 

conditions 


The  constant  K remains  undetermined  and  its  value  must  be  adjusted  for 
correct  temperature  at  prescribed  time.  According  to  the  author,  this 
method  requires  "great  labour". 


For  solidification  in  the  range 
of  the  moving  plane  given  by 


we  assume  a position 


(27) 


where  b,  , and  , are  to  be  determined  from  the  condition  that 
the  temperature  at  the  moving  plane  is  equal  to  the  melting  point  of  the 
medium.  For  continuity  of  the  two  functions  of  *T  , (11)  and  (27) 
and  their  first  derivatives  at  x = x..we  make 


(28) 
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We  determine  b and  from  the  followings 


^ • -fc+ * 4. 

I * ”i»  +•  2 

The  second  equation  of  (<.0)  refers  to  conplete  solidification  at 

J ^ ®T*  X U a J n A vu4  <Sm  M + « +U/-»  1 1 /M.r!  Mr*  A /-\v»  f n ^ 


(29) 


a definite  *Z^  , to  be  obtained  according  to  the  following  considerations. 

While  the  source  of  heat  moves  from  xzO  to  x ~ x# 
temperature  (Vj_  );  is  obtained  from  (18),  and  is  a constant  at  any 
point  in  this  range.  At  the  same  plane  source,  cn-V  will  vary  from 
zero  to  a known  value.  Solidification  is  cor.plete  when  the  source 
movin':  from  x = 0,  and  the  source  moving  *>om  x = 2a  > merge 
at  x = a . At  this  position  we  will  obtain  *“ 
regardlessT  of  their  value  at  any  other  point. 


( VL  )-„•  ( T-l  )t 


As  the  plane  source  moves  from  x 


( V l ) may  be : 


- x_  to 

0 


x = a , temperature 


(a) 

(b) 

(c) 


conti nuously  inc  reasin  .g 
continuously  decreasing 
constant  * WV7t 


-o 


4/a«.)  (i-eUj) 


The  rate  of  solidification  for  (c)  can  be  easily  obtained,  as  the 
temperature,  due  to  the  original  supply  of  heat  at  x - a and 
T » T-  will  be  given  by 

-f2n4l)s'Tt*T<,/4 


1 

Zn+I 


fcn+ll 


E 

” ^ n»c 

= v - 4' e*/ e [i-e (4jJ 


(30) 


The  value 
numerically. 


of  (30)  may  be  found  graphically  and  checked 


The  functions  of  T , corresponding  to  (a),  (b),  and  (c),  will 
be  obtained  by  substituting  in  the  second  equations  of  (29)  a value 
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(a)  greater  than  in  (30) 

(b)  smaller  than  in  ( 30) 

(c)  as  in  (30) 


The  speed  of  the  moving  planes  for  (a),  (b),  and  (c),  will  be 

1 br  «lx  _ 4*  X.< 

T.iTiT, 


(31) 


and  for  each  case  will  depend  on  and  L ^ 

In  any  numerical  example,  the  lowest  rate  of  solidification  will 
be  obtained  for  case  (a)  and  the  highest  rate  for  case  (b). 


The  rate  of  solidification  is  given  also  by  the  rate  of  conduction 
of  heat  at  the  source  * , or 


Kf>. 


(d&), 

14*  solid 


molten ) 


(32) 


* Ll/tfitfoot,  loc.  cit. 
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(Vl) 


Cag.i  (b),  ranr* j <rtod  in  Fi<\5.  shows  th<?t,  if  (V^)f  decreases, 

the  -raDhical  rate  will  be  lower  than  the  computed  rate.  Their  difference 


tends  to  zero  when  ( V ^ 


-►(vL  ) 

1 


X ■ X, 


♦Dotted  lines  on  Fit  k,  Z,  and  6 are  temperature  distributions  (V^_)  . 

Solid  lines  are  drawn  tangent  to  distributions  aw  the  sources. 
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, I 


Fi-  5 

In  cas  : (c),  figure  6,  graphical  and  co-rut  jd  rates  v:.l  11  coincide, 
ns  Lae  -ovin,  I'lanos  x and  2a  - x an > roach  x - a , the  rate  of 
solidification  decreasef and  lends  to  a constant.  i constant  temperature 
(V  ^ )(  , at  .in;-  position  of  the  nlnoo  source  -ovin  r from  x = 0 (or 
( at  cii  • ■dan-;  source  r.cvinp  Cron  x ~ _Pu)  is  the  correct 

a ' -'I- , j on. 


Ve  may  s"  ate  tvrt  in  a freezing  mad? um.  'nitially  a'  temperature 
^ , and  with  its  bourn  iny  surface*-  x » 0 , and  3t  : 2a  , kept  at  zero,  the 

positions  of  the  surface;,  of  3 ana  ration  of  sona  and  liquid  pnases  are 
given  by  the  two  functions  of  *C 


* 2 4.  VT  o«7(T, 

£ . - b + 2 ^ T.  £ 

The  constants  and  the  ranges  of  T of  (33)  will  vary  for  each 
numerical  case  and  are  to  be  determined  from  the  data  of  the  problem 
according  to  the  conditions  given  in  this  section. 


In  the  forms  of  (33),  tha  thermal  di ffusivity  k,  included  in  the 
dimensionless  variable  T*  ^ T > not-  undetermined.  ..e  assume 
k constant  and  equal  to  the  average  di  ffusivity  of  the  medium  in  the 
temperature  ran  ye  V to  ^ , i.e.  in  the  molten  material.  This  constant 
is  to  be  dei  ’rmi  • ;d  ext « riraent ally,  and  methods  o f measurement  may  be 
cased  on  solutions  given  in  the  following  sections. 


The  complete  exact  solution  of  the  problem,  for  the  special  case 
9 - V (i.e.  initial  temperature  is  at  melting  point  of  the  medium) 
is  given  by  Lightfoot.  In  this  case  solidification  proceeds  from  each 
boundary  £ — 0 and  at-  = 2a.  , independently,  and  trie  position  of 
the  moving  plane  x^  (or  2a  ~ x) . until  solidification  is  complete,  is 
given  by  „ 

i 1 2 

« 

(34) 

.here  Jf  is  obtained  from  (2)  as  for  a semi-infinite  medium. 


In  any  numerical  example 


^•v 


, we  can  write 


(35) 


and  we  verify  that  for 


as  in  (35)  we  will  obtain 
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©o 

2 


sin 


Jin±lhL 

2 


sV-Ak3®  «.*V<)  (■-•«)) 


(36) 


Some  numerical  solutions  for  media  having  the  sane  melting  point  are 
shown  in  Figure  7.  The  initial  temperatures  f and  ratios  » are  ad- 
justed for  equal  value  of  constan 


The  Sphere  and  the  Infinite  Cylinder 

In  the  solution  of  the  problem  of  solidification  in  the  finite  region 
2a  , the  planes  of  separation  of  solid  and  liquid  phases  have  been  regarded 
as  moving  plane  sources  of  heat.  Their  speed  was  obtained  from  the 
condition  that  the  temperature  at  the  sources  is  equal  to  the  melting  point 
of  the  medium.  The  temperature  due  to  the  original  supply  of  heat  is 
evaluated  from  a known  solution  while  the  temperature  due  to  latent  heat 
of  fusion  is  determined  from  analysis.  Solification  in  a sphere,  or  in 
a cylinder,  gives  rise  to  spherical,  or  cylindrical  surface  sources  of 
heat,  and  the  forms  of  , for  these  cases,  are  similar  to  (24)  or  (25). 

As  for  moving  Plane  surfaces,  correct  evaluation  of  the  series  of 
integrals  can  not  be  obtained  by  available  methods.  The  temperature 
(VL  ),  at  the  surfaces  of  separation  of  solid  and  liquid  phases  will  be 
determined  from  analysis. 

The  known  solutions  VW  in  the  slab  of  thickness  2a,  in  the  infinite 
cylinder  £ > jr»  ^ 0 , ana  in  the  sphere  4 > O , “Initially  at  f 

and  with  boundary  surfaces  kspt  at  zero,  expressed  in  terns  of  the  dimension- 
less variable  T * *%  T > are  as  follows: 
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M y* 

Tt  ~ 


n*o 


I 

2 n+t 


sin 


fcn+lbTX 

z* 


?ir  JcCraal  e 

•nJ.f-oO 


-ten*!)1! Xzfy4 
e 


24 ^ (-l) 

**r  £ro  n + l 


sin 


/n+l)nr 

4 


-(n+/)zn% 

e 


(37) 


(38) 


(39) 


where  the  subscripts  1,  2,  and  3,  refer  to  the  slab,  the  cylinder  and 
the  sphere,  in  that  order. 


The  temperature  distributions  obtained  with  (37),  (38),  and  (39), 
for  the  same  medium,  with  dimension  a and  initial  temperature  f , 
the  same  for  the  three  regions,  will  never  result  as  in  Figure  8,  regard- 
less of  Tj  , TV  , and  Xi  . Identical  distributions  in  any 
miven  range  X ■ 4 mJTt  *4-X-  and  6<X^4  , may  be  obtained 

for  particular  values  of  , and  tj  , as  in  Figure  9 . 


Ptfl 


r *0 


Fig,  9 
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In  the  problem  of  solidifi cation  we  assume  that,  in  all  three 
regions,  the  surfaces  of  separation  of  solid  and  liquid  phases  Are  at 
equal  distance  from  the  boundary  (we  may  have  ^ T1  ) and  for 

these  surfaces  we  write  9 

v,  t vu « v 


(40) 

At  these  surfaces  we  obtain  also 

(v*),  =(v,)l=(vf,)3 
<vL),  - (VJ* » (vL>, 

(41) 

where  ( V^)f  , in  thi3  instance,  is  the  temperature  caused  by  all  tue 
heat  of  fusion  flowing  in  the  slab  2a. 

We  will  prove  briefly  that  (41)  is  a necessary  condition  in  the 
solution  of  these  problems. 

We  consider  the  slab  and  the  infinite  cylinder,  soli:  at  unit 
distance  from  the  boundary,  or  at 

x » a.-p  «l 


(42) 

and  at  the  surfaces  of  separation,  we  make 

(V,),>(V„)t  (43) 


The  inequality  of  (43)  means  that,  in  the  solid  part  of  the  slab, 
the  temperature  distribution  ( V#  ),  would  be  higher,  at  all  points,  than 
the  temperature  distributions  ( j^in  the  solid  part  of  the  cylinder. 
During  solidification  of  unit  volumes  of  medium  adjacent  to  the  external 
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surfaces,  the  total  flow  of  original  heat*  per  unit  area,  would  be 
greater  at  the  cylindrical  surface  41  >1"  j|fc  - I than  at  the 
corresponding  plane  surface  K*  Mml“  . During  solidification,  in  any 
region,  original  heat  and  latent  heat  of  fusion  flow  to  the  boundary 
concurrently  and  are  subject  to  the  same  boundary  conditions.  With  a 
greater  total  flow  of  original  heat  we  would  obtain  also  a greater  total 
flow  of  latent  heat  of  fusion;  the  temperature  distribution  (V^  )#  in  the 
solid  part  of  the  slab  would  be  higher,  at  all  points,  than  the  temperature 
distribution  ( ^ in  the  cylinder. 

Assumption  (43)  leads  to  the  conclusion  that,  either  in  the  slab  or 
the  cylinder,  the  temperature  at  the  surface  of  separitiop  of  solid  and 

liquid  phases  would  result 

V,  +VL  * V 

(w) 

The  same  conclusion  will  be  reached  reversing  the  inequality  in  (44) 
and  for  subsequent  elements  of  volume,  until  solidification  is  complete. 

The  same  result  will  be  reached  comparing  solidification  in  a slab  and  in 
a sphere . 


This  conclusion  opens  the  wav  to  the  so’ution  of  the  problems  of 
solidification  previously  discussed  and  appears  to  have  been  overlooked 
in  past  work  on  the  subject. 

Since  the  solution  of  the  problems  of  solidification  for  the  slab, 
the  infinite  cylinder,  and  the  spher*,  must  satisfy  condition  (41),  the 
solutions  for  the  cylinder  and  the  sphere  can  be  derived  from  the  solution 
obtained  for  the  slab. 

For  solidification  in  the  slab,  at  a given  distance  */*  , we 
evaluate  % from  the  two  forms  of  (33)  and  solve  (37)  for^rhese  values  of 
T and  In  (38)  and  (39)  we  substitute  y t and 

as  obtained  from  (37).  The  values  and  Tj  given  by  the  (38)  and  (39) 
can  be  obtained  graphically  and  checked  numerically. 

ESgunlO  shows  thet,  initially,  the  rates  of  solidification  in  all 
three  regions  coincide.  After  initial  periods,  the  rates  of  solidification 
of  a cylinder  and  a sphere  will  assume  constant  values;  in  a sphere  more 
rapidly  than  in  a cylinder. 


5,  The  Slab,  the  Infinite  Cylinder,  and  the  Sphere,  Initial  temperature 
constant  and  radiation  at  the  surface  into  a medium  at  constant 
temperature 

A.3  in  the  previous  cases,  the  solutions  of  these  problems  are  derived 
from  the  solution  obtained  for  solidification  in  a slab,  with  initial 
temperature  constant  and  bounding  surfaces  kept  at  zero. 

In  a medium,  initially  at  ^ and  with  radiation  at  x and 

x = +a  into  a medium  at  zero,  the  temperature  distribution  Va  is 
given  by 


tS„+A  + A*]co»A 


(45)* 


where  A — ah,  h being  the  heat  transfer  coefficient  of  the  surface, 
and  A>  At-~An  the  positive  roots  of 

(46) 


The  temperature  distribution  obtained  with  (46)  for  any  value  of  n and 
*U  and  the  temperature  distribution  obtained  from  (37)  for  the  slab, 
for  any  value  of  T , will  never  cross  each  otter,  as  in  Figure  8.  The 
difference  between  temperature  gradients  at  corresponding  points,  through 
the  entire  region,  will  always  be  ^ 0,  or  always  £0  . We  conclude  that, 
when  in  the  region  • a < x < +a  , at  constant  initial  temperature  and 
radiating  at  x«-a  and  at  x*  + a into  a medium  at  cons tant  temperature, 
we  find  the  surfaces  of  separation  of  solid  and  liquid  phases  at  the 
same  distance  from  the  boundary  as  in  the  slab  of  (37),  condition  (41) 
will  be  satisfied.  The  solution  of  the  problem  for  this  case  is  obtained 
by  the  method  outlined  in  section  4. 

Condition  (41)  applies  also  to  an  infinite  cylinder  and  to  a sphere, 
initially  at  constant  temperature  and  cooled  by  radiation  into  a medium 
at  constant  temperature.  The  two  problems  can  be  solved  by  the  same 
method  and  forms  of  VA  for  these  cases  can  be  found  in  textbooks. 


* Carslaw  and  Jaeger,  Conduction  of  Heat  In  Solids,  p.  100 
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The  Finite  Cylinder  and  the  Parallelepiped 

In  the  cylinder  of  Figure  11  initially  at  constant  temperature  ^ and 
with  boundary  surfaces  kept  at  <,ero,  the  surface  of  separation  of  solid 
and  liquid  nhases  will  be  at  a point 


Fig  11 

when  (V^Lis  equal  to  the  temperature  y*  obtained  in  the  slab  of 

(37)  at  thtf  surface  of  separation  » jf  .In  general, 
solidification  in  a finite  cylinder  occurs  with  the  hi  gher  of  the  two 
temperatures  Vd  obtained  for  the  slab  at  the  two  positions  given  by  the 
coordinates  of  the  point  as  in  (47). 

The  known  solution  Vd  in  a finite  cylinder,  initially  at  (fr 
-nd  with  boundary  surfacesrkept  at  zero,  expressed  in  terms  of  the 
dimensionless  variable  , is  as  follows: 
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< cj  f-.-’ud- 


? J mJC 

*-»  Zm+i 
m*o 


tea  lZa±!h*^(tmH) 
zl 


(48) 

Solution  for  a finita  cylinder,  initially  at  ^ and  coo lad  by 
radiation  into  armedium  at  constant  temperature  is  given  in  textbooks. 

Solidification  in  parallelepiped  is  solved  by  the  same  method. 

Fiiyire  32 shows  surfaces  of  separation  in  a finite  cylinder.  The 
numerical  data  are  the  same  as  for  (b)  oi  Figure  ?. 


Iuclosures : 

' - Figure  2 

. Figure  7 
. Figure  10 
« Figure  12 
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I bbbbE  Ebbbb  bbbbb  bbbbb ibbbbbbbbb bbbbbbbibb bbi 

■BBBBBBBBBB BBBBB BBBBB BBBBB BBiBB  BBBBB BBBgB  BBI 
IbbbbEbbbbbbbbbbbbbbbbbbbbbbbbbbbi 


^HBBBBBBBBBBBBBBg  BBBBBBBHBM 
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IBBBBBBBBB BBBBB BBBBB BBBBBBBBBB BBBBB BBBBB BBBBB BBBBB BBBBBBBBBB 
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